Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 1 February 2008 (MN WT^ style file v2.2) 



Waves and instabilities in dissipative rotating superfluid 
neutron stars 

T. Sidery^, N. Andersson^ and G.L. Corner^ 

^ School of Mathematics, University of Southampton, Southampton S017 IBJ, United Kingdom 

^ Department of Physics & Center for Fluids at All Scales, Saint Louis University, St. Louis, MO, 63156-0907, USA 
1 February 2008 

ABSTRACT 

We discuss wave propagation in rotating superfluid neutron star cores, taking into 
account the vortex mediated mutual friction force. For models where the two fluids 
co-rotate in the unperturbed state, our analysis clarifies the role of chemical coupling 
and entrainment for sound and inertial waves. We also investigate the mutual friction 
damping, providing results that demonstrate the well-known fact that sound waves 
propagating along a vortex array are undamped. We show that the same is not true for 
inertial waves, which are damped by the mutual friction regardless of the propagation 
direction. We then include the vortex tension, which arises due to local vortex cur- 
vature. Focussing on purely transverse inertial waves, we derive the small correction 
that the tension induces in the wave frequency. Finally, we allow for a relative linear 
flow in the background (along the rotation axis). In this case we show how the mutual 
friction coupling may induce a dynamical instability in the inertial waves. We discuss 
the critical flow required for the instability to be present, its physical interpretation 
and the possible relevance it may have for neutron star physics. 

1 INTRODUCTION 

Our understanding of neutron star interiors currently relies almost entirely on observations of the "surface" . In some cases we 
have upper limits on the temperature, which can be combined with an estimated age to yield information about neutron star 
cooling. This in turn depends on the interior physics, e.g. whether the core contains superfluid components or not (Page et al 
2004). Other evidence comes from the way that the neutron star interacts with its environment, e.g. how the magnetosphere 
affects the spin-down of an isolated radio pulsar. It is, however, much harder to link this information to the properties of the 
interior. It is also difficult to draw definite conclusions about the nature of the core fluid from bulk quantities such as mass 
and radius. 

The currently most potent tests of our theoretical ideas are provided by observed crust oscillations in the tails of magnetar 
giant flares (Strohmayer & Watts 2006), and glitches in the spin-down of radio pulsars (Lyne, Shemar & Smith 2000). It seems 
plausible that the crust motion in the magnetar events will to some extent depend on the core physics, e.g. whether the magnetic 
field penetrates the core or not (Glampedakis, Samuelsson & Andersson 2006; Levin 2007). Meanwhile, the glitches remain 
the strongest indication that the core contains partially decoupled superfluid components that may (probably following the 
onset of some instability) transfer angular momentum to the crust. In order to improve our models of these events we need to 
understand the dynamics of large-scale superfluid systems. One key question that must be addressed if we want to be able to 
compare our models to real data concerns how energy is dissipated in the system. Consider the (relatively simple!) case when 
the internal fluid is a mixture of superconducting protons, electrons, and superfluid neutrons. The equations of motion for 
such multi-constituent fluids have been formulated (Prix 2004) and constrained to a three fluid model (Andersson & Comer 
2006). The three fluids are neutrons, entropy and charged particles. It has recently been argued that for this system there are 
19, more or less unknown, dissipation coefficients (Andersson & Comer 2006). This problem is clearly much more intricate 
than the standard single fiuid case, where one need only worry about shear- and bulk viscosity. 

A natural way to gain insight into the nature of a fluid system is to carry out a local analysis of wave propagation. 
Such a plane-wave study should provide a better understanding of energy dissipation in the system, and could perhaps also 
help constrain the different parameters. Given our current understanding it is natural to divide this effort into a number of 
steps. This flrst study is focussed on the two-constituent model which applies in the low-temperature limit. In this limit, we 
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know from results for superfluid Helium that the main dissipation mechanism is due to the presence of rotational vortices. 
The vortices induce a mutual friction between the two constituents. We have recently discussed the associated force for 
neutron stars (Andcrsson, Sidery & Comer 2006), including the important effect of cntrainment. The entrainment is due to 
the multi-fluid aspects of the problem. When the energy of the system is allowed to depend upon the relative velocities of 
the constituents, each constituent momentum is modified in such a way that it is no longer parallel with the corresponding 
transport velocity. The model has since been extended to include curved vortices and a suggested form for the average force 
in a polarized state of turbulence (Andersson, Sidery & Comer 2007). 

Since we now understand some of the effects that should be present in a dissipative superfiuid neutron star, it is interesting 
to ask how they affect motions in the fluid. Hence, we will study local oscillations in the fluid in various situations. As the 
problem soon becomes intractable unless one introduces a number of simplifying assumptions, we aim to simplify the problem 
to the point where we can carry out the analysis analytically. Nevertheless, we believe that our results shed new light on 
the dynamics of these systems. Through various dispersion relations and the properties of the associated waves, we get a 
better understanding of the role of entrainment and the nature of two-fluid inertial waves. We also learn how mutual friction 
damps both acoustic and inertial waves. Finally, we account for the vortex tension and demonstrate how the mutual friction 
may induce an instability in the inertial waves. This is the first demonstration of the so-called Donnelly-Glaberson instability 
(Glaberson, Johnson & Ostermeier 1974) in a neutron star. We show that the instability belongs to the general class of 
two-stream instabilities (Andersson, Comer & Prix 2004) , and discuss how it may trigger turbulence in the neutron star core. 



2 THE SINGLE FLUID PROBLEM 

Before considering the superfiuid problem we will discuss the simpler case of a one fluid system. This is a fairly standard 

analysis, but it provides a number of results that will be of use later. Most importantly, it allows us to discuss the nature of 
inertial waves. We also have an opportunity to establish the notation and the approach to the plane-wave problem that will 
be used throughout the paper. 



2.1 The plane- wave approach 

The equations of motion for a single fluid in a rotating frame are, in a coordinate basis, 

{§i '"'^') '^^'i''^' + ^'('^ + *) = . (1) 

The velocity is given by Vi, Q.i is the angular velocity of the frame and jl is the chemical potential of the fluid per unit mass. 

The gravitational potential is represented by Throughout our analysis we will assume that the effects of gravity can be 
neglected. In effect, we make the Cowling approximation by ignoring variations in $. The continuity equation is 

^^+V,{pv')=Q, (2) 

where p is the mass density. Assuming that the background motion is represented by solid bulk rotation, the corresponding 
velocity field vanishes in the rotating frame. Perturbing the system we then have 

5v^ = t;ie'(-'+'=.-') , (3) 

i.e. the perturbation is written in the form of a plane wave, with frequency u! and wave vector ki. Similarly we have for the 
perturbed density 

5p = pe'('"*+'=J'"') , (4) 

The amplitudes Vi and p will be taken to be constant throughout the analysis. With these assumptions, the continuity equation 
becomes 

ujp + pkjV^ = , (5) 

where p now represents the background density. 

To perturb the equations of motion we assume a one-parameter equation of state. Representing the equation of state by 
an energy functional E = E{p) we then have 

dE djl 2I0 

''=d^^^''=d^^P = ''^p^P' 

where the sound speed is defined as 
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It now follows that 

iuJVi + 2eijkQPv'' + ikiC^-p — Q . (8) 
Substituting the continuity equation (5) into (8) we finally find that 

iujVi + 2eijk^^v' — ic^ki = . (9) 

UJ 

2.2 Finding the dispersion relation 

The derivation of the required dispersion relation is straightforward. It is preferable to work with scalar equations. Hence we 

first contract (9) with fc*. Rearranging the result we have 

i{k'vi) (^uj - ^ j - 2n'eijkk'v'' = . (10) 
Contracting (9) with Q.ie^^'^kj then leads to 

iLj{Q*eijkk'v'') + 2Q''{kjv') - 2{kjQ'){QiV*) = . (11) 
Finally, contracting (9) with gives 



■jQ' = ^{kjn'){vik') . (12) 



Defining the angle 9 such that kjil^ = kflcos6 we can now use (11) and (12) to substitute for tijk^^k^v'' and fljv^ in (10). 
Thus we get 



{v,kn 



+ 4 — 1 — cos " 



. (13) 



U) J UJ \ UJ 

Provided that the wave is not purely transverse, in which case we would have Vjk^ = 0, we arrive at the dispersion relation 

- uj"^ [clk^ + 4n^] + 40^ cos^ Ok^cl = . (14) 
In the limit of slow rotation, this quartic in a» has approximate solutions, 

UJ w ±Csk \1 + -^^sm^ 0\ w ±Csk , (15) 

UJ w ±20 cos 61. (16) 

These are the well known results for sound waves and inertial waves, respectively. 

Before we proceed, we need to consider whether the system may admit purely transverse waves. This is important since 
(14) does not apply when Vj^ = 0. In this case, we first of all see from (12) that we must also have parallel with the wave 
vector A:'. The corresponding dispersion relation follows easily by taking the cross product of the Euler equation (9) with f2'. 

This leads to 



40' 



2 - 

uj-^]vi=Q . (17) 



Hence we have w = ±2f2. In other words, we obtain the 9^0 limit for the inertial waves. Thus we learn that, in this limit 
the inertial waves are purely transverse. By returning to the perturbed Euler equations, and representing the solution in a 
Cartesian coordinate system where the z-asas is aligned with the wave vector, we see that we must have Vx = Hvy. This 
solution obviously satisfies VjV^ = 0. 



3 THE TWO-FLUID PROBLEM 

Let us now extend the plane-wave analysis to the two fluid problem. In principle, we expect a doubling of the number of 

solutions. In addition, we are interested in the new features that become relevant when we are dealing with a multi-fluid 
situation. Our focus will be on the entrainment, which represents a non-dissipative coupling between the fluids and the 
mutual friction, which represents dissipation due to electrons scattering off of the neutron vortices. 

As we want to account for vortex effects it is natural to consider the equations of motion in a rotating frame. We will 
assume that the background equilibrium is such that the two fluids rotate together. This situation, which would correspond to 
the two fluids being in chemical equilibrium, is slightly simplifled which keeps the analysis manageable. Having said that, it is 
worth emphasising that we also need to understand what happens when there is a velocity difference in the background. Since 
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relative rotation is required for the standard explanation of the pulsar glitches, one would in principle expect the two fluids in 
a neutron star core to rotate at different rates. If we are primarily interested in dynamics on a short timescale compared to the 
time it takes to re-establish chemical equilibrium we have the freedom to consider different background velocities. This general 
problem is, however, quite complicated and it makes sense to first consider the case where the two fiuids rotate together. 

As discussed by, for example Prix (2004); Andersson & Comer (2006), the required equations of motion can be written 
(in a frame rotating uniformly with angular velocity f2') 

{§i + ''^'^^■) + ^-""^'^^ + Sxwf Vit;i + Vi($ + Ax) + 2eijkfi'v^ = /f . (18) 

The constituent indices x and y (x y) label the fluids, and will later take the values n and p. The former represents the 
superfluid neutrons and the latter a charge neutral conglomerate of protons and electrons. The relative velocity is denoted 
by w^^ = Vi — . The force /f on the right-hand side represents a dynamical coupling between the two fluids. We will focus 
on the case when it arises due to electron scattering off of the magnetic fields associated with the rotational vortices in the 
neutron superfluid. Then we have, assuming that the vortex array is straight (Andersson, Sidery & Comer 2006), 

fi = —UvB' eijkK^ w^y + —nvBeijkK^ e'''"' Kiw^ . (19) 

Px Px 

In this expression n„ is the vortex number density per unit area, and at the macroscopic level we have (Andersson, Sidery & 
Comer 2007) 

rivKi = eijfcV^(Vn + enWpn) ■ (20) 

This means that for a straight vortex array, representing bulk rotation, we would have 

rivKi = 2n° + 2en(nP - n^) . (21) 

In our chosen background conflguration the two fluids rotate together, so the second term in the above expression vanishes. 
As usual, mass conservation requires that 

^+V,(px«i)=0. (22) 

We consider small perturbations away from a background where the two fluids are at rest in the rotating frame. That is, 
we use 

^+V^p^Sv])=0. (23) 
Making also the Cowling approximation (5$ = 0), the perturbed equations of motion become 

^ {5vt + exSw^) + ViSpx + 2eijk^'Sv^^ = 5ff , (24) 
where 

Sfi = —rivB'eijkK^ Sw^y + —nvBeukk^ e'''"^ kiSw^ . (25) 

Px Px 

It should be noted that setting the background rotation rates equal greatly reduces the complexity of the perturbed mutual 
friction force. 

We next assume that the perturbations can be represented by plane waves. Then we have 

5px = Pxe'('^*+'=^=^'' , (26) 

where vf and px are assumed to be constant (that is, they vary on lengthscales much longer than the wavelength). The 
continuity equations then become 

kjvi = -^uj. (27) 

Px 

To make progress we need a representation of the equation of state. In our formalism (Prix 2004; Andersson & Comer 
2006), we need to provide an energy functional E from which the chemical potentials follow according to 

dE 



Op: 



(28) 

The entrainment Ex is similarly determined as 



2a ^ dE 
Sx = — , where a 



Px ' dw^ 



(29) 
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The basic idea is that if the system is isotropic, then the energy function can be constructed from the various scalars that can 
be formed from the dynamical variables px and w^- This means that one would generally expect to have 



d/Ux = — opx + — opy + „ „ ow,, 



(30) 



In the present case, where the two fluids move together in the background, the last term will vanish since 
Hence, we can use 

5/ix = ^—Spx + ^—Spy = P-kkSPk + Axy^Py . (31) 

OPy 

Later, when we allow for relative flow in the background, this form for the perturbed chemical potentials will no longer be 
generally valid. We will nevertheless use it in order to simplify the analysis. In practice, this amounts to assuming that the 
equation of state is "separable" in the sense that it can be written^ 

E = f{pn,pp) + .f;(wVx) • 



Under the above assumptions, the momentum equations become 



iu){vl + Exwf'') - iki{ji^y,p^ + Axypy) + 2eijk0.^v^ = 2— B'eijkii' w^y + 2—Beijkk^ e'''"^Qiw^^ 

Px Px 

After using the perturbed continuity equation (27) and rearranging we arrive at 
iu}{l - s^)Sr + 2eir^' - 2^B'eij"'Q' - 2^Beijkk'e''^"'Qi - ikiji^ 



(32) 
(33) 



pKL Px " ' U> 

^ I o P^ 12' ^r\j I o P^ K? "j klmr) -1 ~ Py J ■ 

luie^Oi +1 — B eij i2 +2 — Beijki^^ ili ikifi^y — k 
Px Px ijJ 

Let us now introduce the speeds of sound as (Andersson & Comer 2001) 

2 _ 
Cx — Px/^xx , 

and represent the "chemical coupling" by 

Cx — Px/-ixy — PxPyx • 

Since the partial derivatives commute we have 
C - 

i_-p — i^n ■ 
Pn 

Then we have 



. 



(34) 

(35) 
(36) 
(37) 



tuj{l — SxjOi + 2eij il-^ — 2 — B tij il^ — 2 — Btij^K' t Vli — iki — k 



px ' Px 

Px " px ' W 

The dispersion relation for wave propagation is encoded in this equation. As in the single fluid case, it is preferable to work 
with scalar equations. Hence, we contract (38) with fc,. This gives 



iwExC + 2^B'eij"'^l' + 2^Beijkk'e''^"'(li - ih^k"" 



. 



(38) 



iw(l - Ex)*;"" + ( 2 - 2^B' ] k'eir^' - 2k' ^Btijkk' e^'^^Q-i - ik'^k 



„2 
,2Cx , 



Px 

iwExfc™ + 2^B'k'eij"'n' + 2^Bk\^jkk' e'^'^^ni - ik'^k" 

Px Px 



px 



2 Cy 



(39) 



We obtain a second scalar equation by contracting (38) with Q.it'^ kj 



Px 



ia;(l - e^)Q.ie"'"'kk + I 2 - 2f^B' ) Q.it"''kkeij"'Q.' - 2f^Bnie""-k^trjkk' e'^^'^'Q.i 



Px 



lioeAe'^'^kj +2^B'flie''''kkei,'^Q^ +2^Bni€'''''kgerjkK'e''^^^^ 

Px Px 



A third equation follows from the contraction of (38) with 



= . 



(40) 



(41) 



^ It is worth noting that if one wanted to propose a "realistic" separable equation of state, then one would have to think more carefully 
about the units of the two terms. In other words, our proposed form for E should be seen as a mathematical construction motivated by 
the fact that it simplifies the analysis. 
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Although in principle straightforward, the formulation of the dispersion relation is still rather messy. To facilitate its 
construction it is useful to introduce some further notation. To motivate this, let us write out the two equations (41) explicitly. 
We have 



If we define 

J. / 1 - En 



En 
1 — En 



then (41) can be rewritten as 



M = 



iivepn"" - ikiil'—k" 



c ^ 



= , 
= , 















Vk = 1 





(42) 
(43) 

(44) 



ki^ 



■^vn = jivk ■ 



As a note for future reference, we will also use 



l-£. 



1 - En 



^ = det ^ = 1 — Sn — Ep 



Q.'-tijkk^vl 

^ ^ijkk'^ '^p 



We can now rewrite (39) and (40) as 



and 



Pp 



jl + 2 — fipo \ Vk - [21 - 2 — po ) V 



Pp 



,BQ. 



pp 



,B' 



Pp 



iuj$, + 2 po + 2/ - 2— p„ n'vk - 27 - 2— po {kjQ')vn = 



2-^{kjQ')poVn = 



,B' 



pp 



Assuming that 1 — En — £p 7^ 0, and substituting (45) into (47) and (48) we get 



B 



iuj^-i — p, + 2 — -po n - 

OJ ppil ' 



Vk 



B' 

21 - 2— Po ]v, = , 
Pp 



and 



We now define 



B' 

21 - 2— Po 
pp 



Vk ■ 



zero 

Po 



Pn Pp 
Pn Pp 



P = 



Pp -pp 
— Pn Pn 



such that p1^''° Pq — 0. The matrix on the left hand side of equation (50) then has the inverse. 



-uj^i + 2iujQ.B 1 + 



pn 



Pp 



Using this in (49) we obtain 

lUJt, - I — M + 2 ;=rPo 

w ppf2 " \ 



.i!:*^ — 1 I „zero 
lLO(,i +2 Po 

Pp 



B' 

-[21-2— p^ 



) ^zero 



Vk 



iujiC' + 2^pl^ 



Pp J ^2^ _ 2iuQ.B {\ + 
By using PoP^''° = this expands to give 

2 (fc.fi^T 



2/-2^Po ) ( Q. 



2 (fcj^^) ^-1- 
,,,2 ^ ^ 



Vk = Q ■ 



Vkluo^-i-^ + 2—p,[Q 



ui^^ - 2iuinB 1 + ^ 



(4iw$4 l-4^a;^ [pj,^ ipj + ja;d— PoC Po+8— 



V Pp 



Pp 



Po 



4iw^ ^1 ^/I - 4^^ — (po^ Vo)^ V + 

Pp 



26' 
Pp 



Po^ PoS P + 8^Po ^ P 



(45) 



(46) 



(47) 



(48) 



(49) 



(50) 



(51) 



(52) 



(53) 



= . 



(54) 
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In order for us to have Vk 7^ the determinant of the matrix in the curly brackets must vanish. This condition provides the 
dispersion relation for waves in the two-fluid system. 

As in the single-fluid problem, our analysis docs not apply to waves that arc purely transverse. Such waves arc, however, 
not very likely unless we align the wave vector with the rotation axis. In the general case, we see from (49) and (50) that we 
can have purely transverse waves (for which vk = 0) only if also = or if 

det I ia;4 + —Po J = , 

det (21 - —Po] = . (55) 



These conditions lead to 

= w^(l-en-ep)-iwnB( 1 + — 

V Pp 

= 4 - 2S' ^1 + . (56) 

The second condition is extremely restrictive and so a purely transverse wave is unlikely. In fact, for typical neutron star 
conditions we expect B' <C 1 (Andersson, Sidery & Comer 2006) which suggests that purely transverse waves are not possible 
unless the wave is aligned with the rotation in such a way that Ve = 0. Although somewhat contrived, this particular case is 
interesting and we will discuss it in more detail in Sections 5 and 6. 



4 ILLUSTRATIVE EXAMPLES 

In the previous section we wrote down all the relations we need to derive the general dispersion relation for the two-fluid 
problem. It should be clear that, since the generic dispersion relation is a high order polynomial in u), this problem is quite 
rich. In order to understand the solutions it is useful to consider a sequence of increasingly complex model situations. This 
will give us a feeling for how the various parameters in the model affect the wave propagation. 



4.1 No rotation, coupling, or friction 

It is natural to start with the very simplest case, with the two fluids completely decoupled. This model corresponds to an 
equation of state of form 

E = f(nn)+g{np). (57) 

This leads to /Xxy = and £x = 0. In essence, the two fluids are not coupled either chemically or by entrainment. If we also 
assume that there is no background rotation or friction in the system, i.e. let Q = B = B' = 0, then the dispersion relation 
follows from the determinant 



iujl — i — a 

UJ 



. (58) 
This expands to 



which has the non-trivial solutions 

up' = fc'^Cn and uP — k^Cp . (60) 
Hence, we have the anticipated result that the system only supports sound waves, 

UJ = ±kCn , 

UJ = ±kcp . (61) 
It is also easy to show that these waves are longitudinal, as one would expect. 



4.2 Including entrainment 

We can now begin to investigate how various coupling mechanisms modify these waves. Let us first consider the entrainment. 
Then we need an equation of state that depends on the relative velocity. Thus, we assume that 
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E = /(nn) + g{np) + h{wlp) 



(62) 



This is obviously not the general case, but since we want to be able to analyze the problem analytically it is natural to restrict 
ourselves to this class of separable models. It is straightforward to study more generic situations numerically, but the results 
should not differ qualitatively from the ones we discuss here. 

Still assuming that ^1 = B = B' = /Xxy = we obtain the dispersion relation from 



UJ 



= . 



We expand this to get 



1'" z 



1 £p o '^P 



. 



(63) 



(64) 



The solutions are found from 

- tj- [(1 - en)fc^Cp + (1 - £p)fc^Cn] + fc^CnCp = . 

Solving for we have 

= ^ [(1 - £n)fc^Cp + (1 - ep)fc^C^] ± ^ { [(1 - ^n)Cp + (1 - £p)Cn]'^ " 4^CnCp| 

To make further progress it is useful to assume that the entrainment is a small effect and use Taylor expansion in Sx. We see 
immediately from (64) that to linear order in entrainment the frequencies are given by 



1/2 



(65) 



(66) 



±(1 + -en)feCn , 

±(1 + ^ep)'scp . 



(67) 



This illustrates how the sound waves arc affected by a weak entrainment coupling. At this level there appears to be no 
interaction between the two wave speeds. This would be a higher order effect for this equation of state. 



4.3 Chemical coupling 

Let us consider the other way that the two fluids in a non-rotating system may be coupled. In order to see what effect chemical 
coupling has on the waves, we consider an equation of state of form 



= /(nn,np) 



(68) 



The key difference is that the chemical potential of one fluid can now be affected by the population density of the other 

constituent. This requires us to work out 



, 



which expands to give 

4 2,2/2, 1\ , , A I 1 1 Ppn2 
W - W A; (Cn + Cp) + CnCp 

V Pn 

Taking as small and solving for lu, we get either 



UJ = ±kCn 



1 + 



1 + 



Pp 



2pn fc2c2(c2 - C2) 

Pp Cl 



2pnfc2c2(c2-c2)J 

These are still modifled sound waves associated with each constituent. 



(69) 



(70) 



(71) 



(72) 



4.4 Slow rotation 

We now move on to the case of slow rotation. In addition to the sound waves, we expect to find inortial modes. Since we are 
assuming that the two fluids co-rotate in the background the inertial modes are likely to be degenerate. To keep the problem 
simple, we assume that there is no chemical coupling. In practice, we again let the equation of state be of the form (57). Then 
taking B = B' = = Axy = the dispersion relation follows from 
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lul — I — u TT — Ai — cos Bu, 



. 



As before, 6 is the angle between ki and fij such that fc-Tij = fcOcos^. The determinant expands to give 



1 ^ 2 
1 



— 5" + 4 cos t'c. 



fc2 ^ 4j^2 02fc2 2 

1 -Cp + 4 — — cos 6lcp 

Ct)^ Ct)^ CJ'* 



LO-" U)^ LJ* 

There arc clearly two decoupled cases. The two sets of solutions are found from 

= ui'^-w'^ [k'^cl + 4n'^] + 40^/!c^ cos^ Ocl , 

= - [fc^4 + + cos^ 9cl , 

If we for simplicity assume slow rotation, the solutions to (75) are 

/ 20."^ 2 
ui ±fccn ( 1 + ^j2^ sin 6 

UJ ~ ±2n cos 6 , 

while the solutions to (76) arc 
-^2 



±kCry 



20.^ 2 
1 + TTrrr sm I 



k^cl 



±20 cos 61 . 



(73) 

(74) 

(75) 
(76) 

(77) 
(78) 

(79) 
(80) 



The solutions (77) and (79) represent sound waves with a correction due to the slow rotation. Solutions (78) and (80) are 
the (in this case degenerate) inertial modes. The form of the solutions is exactly £is one would expect from the single fluid 
problem. 



4.5 Mutual friction 

The simple cases we have considered so far provide an insight into the different classes of waves that will be present in the 
rotating two-fluid problem. We now want to develop an understanding of how these waves are affected by the mutual friction. 
To do this, it is natural to assume that the induced damping is weak. In the neutron star case, we also expect to have B' -C B 
which allows us to simplify the problem. Andersson, Sidery & Comer (2006) showed that B' ~ , and the "canonical" value 
for S is ~ 4 X 10~*. Since we are assuming that the mutual friction is weak, it is natural to include it as a perturbation of 
the solutions we found previously. 

In order to be consistent we cannot consider the effect of mutual friction without at the same time accounting for rotation. 
Without rotation there would be no neutron vortices in the background and hence no mutual friction. We therefore consider 
the situation when both B and Q can be assumed small (in a suitable sense). To make the analysis tractable we assume that 
Ex ~ Axy = 0. Stricly speaking, it is not consistent to neglect the entrainment here. It plays a central role in generating the 
mutual friction since it is the entrained protons flowing around a neutron vortex that generates the main component of the 
vortex magnetic field (Andersson, Sidery & Comer 2006). Hence, if we neglect the entrainment then we should not have the 
mutual friction either. Of course, the two contributions have completely different effects on the dynamics. As long as we are 
mainly interested in the dissipation the assumptions we make here should be acceptable. 

The equation that we need to solve can be written 

. ^ ik^ . 2B /^2r {kj^'f - 
lul /LI + — prpo ^ I - o M 



■ * ~ I o zero ~ 

't4(j;/x + 8 Pq h 

Pp 



= . (81) 



Before substituting the appropriate solutions for wo in (81) we note that, when written out in full, we need the determinant 

of a matrix of form 



a Bb 
Be d 



(82) 



This means that, to first order in B, the dispersion relation is either a = or d = 0. Using this fact, and assuming that the 
frequency ui will be replaced by loo + 5lo where wo represents one of the undamped solutions from Section 4.4, we need to solve 
either 

\ iwo 1 + — 1 )cl + 2- - ^ \ ' ci - — 1 - 2— + 2i 1 + ^ 
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= , (83) 



z4wo 1 + — + 8 pn - ^ \' 1-2— 



Z4a;0 IH Cn + 8 PnCn 

Wo/ Pp , 



5w-2iBf2 . 2 „^ 107.202 „„„2 a„2 • (^^) 



za;o(l H ) (1 )Cp + 2-— Q. ^ — 2 U ~ 2 h 2i iH 

wo wo wo S2pp \ ciJo ^ } a;,, \ wo wo \ pp 

Linearising (83) we find that the mutual friction correction to the waves associated with the neutron fluid is given by 

[luo + 4n^] [loo - cos^ ed] 
"w* + k^ulcl + m'^ujl - 12k'^Q? cos2 eel 

Let us first consider the inertial modes, i.e. take Wq = 4f2^ cos^ 6. To first order in Q, this leads to 

6uj = iBQ. [1 + cos^ e] . (86) 
The correction to the sound waves follows by taking ujq = k^c^ + 4D,^ sin^ 9. This leads to 

Jw = iBfl sin^ e . (87) 
To find the other set of solutions we linearise (84) and solve for Suj. This leads to 
p„ [wo' + 4^2] [wo' - fe^ cos^ 6>c^ 



5w = 2iBn- 



'Pp W^ + k'^UJ^cl + 402^^2 _ 12fc2f22 cos2 5)^2 • 

For the inertial waves we again use Wq = 4,0,^ cos^ and find that, 

5oj = iBQ.^ [l + cos^e] . (89) 
Pp 

Finally, using the proton sound wave solution = fc^Cp + 4f2^ sin^ we have 

5Lj = iBQ.f^suYe . (90) 
pp 

To summarize the results, we now have two sets of sound waves that are damped by mutual friction. Their frequencies 
follow from 

w = ± (fccn + sin^ 0) + iBVl sin 6 , (91) 

w = ±(fccp + 2a%in^6') sin^6l . (92) 

pp 

There are also two sets of inertial modes. In the undamped case their frequencies are degenerate, but they become distinct 
when we account for the mutual friction. These solutions are 

uj = ±2n cos e + iBQ. [1 + cos^ 6] , (93) 

uj = ±2Qcose + iBQ— \l + cos^e] . (94) 
Pp 

From these results we learn the following. First of all, (91) and (92) show that there will be no dissipation of sound waves 
that travel along the axis of rotation. This is natural since the sound waves are longitudinal and the mutual friction only 
affects motion orthogonal to the vortex array. It is interesting to contrast this with the result for the inertial modes. Prom 
(93) and (94) we see that these waves are always damped. In fact, the effect of the mutual friction is maximal when the wave 
travels along the vortex array. This result is easy to understand from the discussion of the single-fluid problem in Section 2.2. 
Since the inertial waves generally have a component that is orthogonal to the vortex array it is natural that they experience 
damping due to mutual friction. 



5 INCLUDING VORTEX TENSION 

Up to this point we have implicitly assumed that the vortices can be considered straight. In effect, we have ignored the 
tension that arises because of vortex curvature. This tends to be a relatively small effect, so one would not expect our results 
to change much if we account for it. However, it turns out that the vortex tension is important for the instability that we will 
discuss in the next section. In particular, it determines the critical wavelength at which the instability sets in. Hence, it is 
useful to extend our discussion in such a way that the tension of the neutron vortex array is accounted for. This discussion is 
modelled on Hall's analysis of the corresponding problem in superfluid Helium (Hall 1958). By redoing his calculation within 
our formulation we will show how entrainment affects these modes. For simplicity, we will ignore the mutual friction in this 
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section. The derivation of the tension term is provided in Appendix A. Including the relevant contribution, the equations of 
motion in a rotating frame are (as before, we ignore the gravitational potential) 

+ U^Vj ) (v" + £nWf") + £nW?°ViV^ + Vi/in + 2eijfef2''Vn = PUv V j fti , (95) 



dt 



As discussed in Appendix A, we have [recall (20)] 



(96) 



UvKi = —CijkV^pt = €i3fcV^[Vn + EnWpn] , (97) 
THn 



while 

1 — £p 



log - , (98) 



1 — En — £p 1 — En— £p47r \ a() ^ 

see Andersson, Sidery & Comer (2007) for a detailed discussion. 

We have already worked out most of the terms we need to discuss the perturbations of these equations. The only new 

piece is the tension contribution. If wc consider the same background configuration as in the previous sections, then the two 
fluids rotate uniformly at the same rate and wc have n„K' = 2f2' in the background. We then need to work out 

Sft'^'''°'^ = 5[jyn^ti'Vjk,] = u^5{n^K')Vjk, + 2n'V,Sk?^ . (99) 

The first term is easily worked out from (97). The definition also leads to 

dki = ^^(tijkV'6p^-Qid'eimnV"'5p'^) . (100) 
If the background configuration is uniformly rotating, we find that 

Sft~ = -^^'^0 [eu^y'Sp^ - Clid'eimnV^Spl) . (101) 
As in the previous sections we now make the plane wave Ansatz, i.e. we assume that pl^ = Pn exp[i{u)t + kjX^)]. Then 

^Jtonsion ^ r p -fc _ j-^^^^^i^m-ji^^ 1 ^ j-^q2) 

rrin L J 

where we have defined kz = kjO? . 

Since our main interest is to see how the vortex tension affects the various modes that we have discussed previously, it is 
useful to make a further simplification at this point. Wc will concentrate on waves that propagate along the axis of rotation. 
Then kz = \k\ and since fc* is parallel to fi* the last term in (102) vanishes. Hence, we have 

This expression shows that the tension has no effect on longitudinal waves that travel along the rotation axis. In other words, 
the sound waves are unaffected by the inclusion of the tension. The same is not true for the inertial waves. 

Combining the above results with results from the previous sections we arrive at the perturbed equations of motion; 

iu)[Vi + en{vf - Vi)] - ifei/innPn + 2eijk^^ = -Pfc^ejjfcF [Wn + en(Vp - Un)] , (104) 

and 

iu)[vf + ep{vi - vf)] - ikip,pppp + 2eijkfl^Vp = . (105) 

While we could work out the dispersion relation for generic waves in this system, we have chosen not to do this. The 

reason is very simple. As already mentioned, when the wave vector is aligned with the rotation axis, as in the above equations, 
then the sound waves are unaffected by the tension. Given this, it is natural to simplify the analysis by focussing on pure 
transverse inertial waves. For transverse waves we have v^k^ = which leads to px = by virtue of the continuity equations. 
Hence the perturbation equations can be written 

iuiPi = —2Qeijkk^v!^ — Pkltijkk^p^ , (106) 
and 

iujp? = -2^1eijkyvp . (107) 
To derive the dispersion relation we first take the cross product of each equation with fc,. This leads to the relations 
iujsijkypl = 2m^ + vklp^i , (108) 
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and 

iuieijkk^pp = 2f7t;" . (109) 

Recalling the definition of the momenta Px we can solve equations (106) and (107) for eijkk^v^. Inserting the results in 
equations (108) and (109) we have 



[2n + {l-enpk^.f -uj^(l-en) 



= <U> en(l - En) 



and 



1 - 



(1 - £p)t/fc^ 

2n 



1 - £n - 

2 



2-7 2 
£n£pOJ Vk^ 



20. 



+ '"^^ J^^^ - [2n + (1 - e^)9ki]e^uki \ v^^ 



20 

, 2i _,2 



r/l(-i2 /I \2 2i-p 2 2-n 

[4n — (1 — Ep) <^ J^i = 'i' EpWi ■ 

Prom these two relations we see that the required dispersion relation is 

1 - En - 

£n(l - £p)w^ 



(110) 



(111) 



[4^2 _ (1 _ spfu:^] \ [20 + (1 - e^)vklf - w^(l - e„) 



2 2 1 
EnEpW Vk 



20 



20 



= W^Ep < a;^En(l - En) 



1 - 



(1 - Ep)l'fc^ 



20 



20 



- [20 + (1 - e^)Pkl]£^Pkl 



(112) 



In principle, it is straightforward to write down the solutions to this equation. After all, it is just a quadratic in lj^ . Of course, 
the resultant expressions will be so complicated that we learn very little from them. Let us instead focus on two limiting cases. 
First of all, we see that if we neglect the entrainment we have 



[40' - u}'][{20 + vklf - w^] = 
The solutions are obviously 
w = ±(2n + vkl) , 
and 

Lj = ±20 . 



(113) 
(114) 
(115) 



The first solution represents the neutron inertial modes, and the second corresponds to the inertial waves in the proton fluid. 
As one might expect, the former arc affected by the neutron vortex tension while the latter arc not. These modes arc analogous 
to the modes found by Hall (1958) in the case of superfluid Helium. Of course, our calculation adds to the standard analysis 
for Helium by accounting for the entrainment. To get a first idea of how it affects the inertial waves, we can include it as a 
small correction to the above solutions. We then find that, to linear order in the entrainment we have 



a; = ±(2n + i/fcf + 2£nn) , 
and 

uj = ±2(i + £p)n . 



(116) 



(117) 



It should, of course, be emphasized here that there is no physical reason why the entrainment parameters should be small. 
We have simply made this assumption in order to facilitate an analytical calculation. 

To summarize, we have shown how the tension of the neutron vortex array provides a small correction to the inertial 
modes in the neutron fluid. We have demonstrated that this remains true when the entrainment is considered weak and the 
calculation is carried out to linear order. The full solution to the problem, obtained from (112), is likely to exhibit a more 
complex structure. This could easily be investigated via numerical solutions of (112) for some suitable model equation of 
state. At this point we arc, however, not going to discuss this possibility. Instead, we will consider the effect of introducing a 
relative flow on the background configuration. 



6 INSTABILITY OF THE VORTEX ARRAY 

So far we have assumed that the two fluids rotate together in the background configuration. This is a natural assumption 
given that dissipation will tend to damp any relative motion. However, there are a number of situations where one may be 
interested in dynamics that takes place on a timescale shorter than that associated with dissipation. Then one can relax the 
conditions of both chemical and dynamical equilibrium. In particular, one can allow for relative motion in the background 
configuration used in the plane-wave analysis. The question is whether a relative background flow alters the solution we have 
discussed in an interesting way. This turn out to be the case. In fact, by allowing a relative flow we will see that the vortex 



Waves and instabilities in dissipative rotating superfluid neutron stars 13 



array may suffer a dynamical instability. This instability is well-known in the case of Helium, and is often refered to as the 
Donnelly-Glaberson instability (Glaberson, Johnson & Ostermeier 1974). Since the two-fluid model for a superfluid neutron 

star core is completely analogous to the standard formulation for superfluid Helium, it is no surprise that this instability is 
relevant also for neutron stars. In this section we derive the critical relative velocity for this vortex instability, and discuss its 
interpretation. 

In order to keep the analysis tractable we extend the case discussed in the previous section. That is, wc focus on purely 
transverse waves in the case when the wave vector k' is aligned with the rotation O*. In addition, we will assume that it is 
sufficient to consider the dynamics of one of the fluids. In practice, we consider the protons as "clamped" and ignore their 
contribution entirely. This setup is analogous to that discussed by Glaberson, Johnson & Ostermeier (1974) for Helium. In 
that case, the assumption can to some extent be justified since the "normal" fluid is viscous. In our case, this would also be 
true, since our "proton" fluid accounts for the electron component, which will be affected by viscosity. It is not clear, however, 
that the viscous timescale is short enough that the clamping assumption is truly justified. This is an important caveat, but we 
do not believe that relax;ing this assumption would alter our results in a significant way. For simplicity, we have also chosen 
to neglect the entrainment effect here. 

We focus on the perturbed neutron equation in the case when there is a relative flow in the background. To facilitate the 
analysis we first assume that this background fiow is aligned with both the wave vector and the rotation axis. Representing 
the background flow by Vnk^, the perturbation equation can be written 

i{oj + V^kj)vf + 2Qeijkk'v^ = dfi . (118) 

The force on the right-hand side has three contributions. The contribution from the vortex tension remains unchanged from 
the previous section. We also need the mutual friction force. Under the present assumptions, and if we also neglect B' which 

makes sense since it is much smaller than B, we get from (19), 

6fr = Bti^ut"'"- \iCliV:^t'^''k^vl + iCl^ti^^k^vlV:^ + 2n'Clivl\ . (119) 

Finally, we also want to account for the contribution to the mutual friction from the self-induced flow. This is a small 
contribution, but it is natural to included it if we are considering the vortex tension. As long as the background flow is 
uniform, this term can be written, cf. Eq. (36) in Andersson, Sidery & Comer (2007), 

fY"^ = -vBn^eijkK'k^Vik'' . (120) 
Perturbing this we arrive at the contribution 

dfr'' = iyB{Cl'ki)eijkCl^e''"'^kmK ■ (121) 

Putting all this together, we consider an equation of form 

iiuj + Kfc.)v° + 2Q.eijkyvl = -{vkl - iBV^k,)eijkk^ - B{2^1 + iyk1)vf . (122) 

Taking the cross product of this equation with k' and combining the two equations we find that the dispersion relation is 
simply 

+ KA;^ - iB{20, + vki)f = [2fi + vkl =F Kfc^]^ - (123) 
That is, the inertial waves in this system must have frequency 

uj + Kfc^ = ±(2n + vkl) + iB{2Q. + vkl =F V^k^) . (124) 

Given that our assumed time-dependence is exp(iwt) this expression shows that the solution corresponding to the upper sign 
will be exponentially growing {oj has a negative imaginary part) when 

2€l 

V^>—+uk,. (125) 
kz 

In other words, for any given wavevector kz there exists a critical relative flow above which the wave is unstable. Of course, 

we see from (125) that the critical flow must be large both in the limits of large and small kz. If we are interested in the 
critical flow at which the instability first sets in, then we simply need to find the minimum of the function on the right-hand 
side of (125). Thus we need 

Inserting this in the expression for the critical flow we see that the system will have unstable waves when 



y„ > Vc = 2V2n;^ 



(127) 
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This is exactly the condition derived by Glaberson, Johnson & Ostermeier (1974) for the Helium problem. 

Even though we have not attempted the general problem, without assuming that the protons are clamped, we have 
relaxed some of the other assumptions. In particular, one does not have to assume that the waves are purely transverse. The 
interested reader can find a more general discussion in Appendix B. 

Let us now see if we can understand the nature of this instability better. To do this it is helpful to consider the phase- 
velocity of the waves. Recall that in the present problem setup, a constant phase would mean that 

Re ojt + kzZ = constant — » ffp = ~—Tr~ ' (128) 
where ap is the phase-velocity. Hence the phase-velocity of the inertial waves is 

fTj, = K T (2n + vkl) . (129) 

Comparing this to the condition for the critical velocity we immediately see that the instability sets in through the waves 
that propagate in the direction opposite the background flow, (for a suitably small Ki). The critical point is simply identified 
as (jp = 0. The interpretation of this condition is that a wave that is originally seen as travelling downwards (relative to 
Vn) is dragged upwards by the flow and becomes unstable when its direction of propagation changes (according to a fixed 
observer). This condition is typical for a two-stream instability. We have previously considered this class of instabilities for 
neutron stars, see Andersson, Comer & Prix (2004) for a discussion and a list of relevant references to the plasma physics 
literature. A two-stream instability typically requires two identifiable flows and some coupling between them. In our previous 
discussion of such instabilities for neutron stars, we focussed on chemical coupling and the role of cntrainmcnt. We now see 
that the instability can also be caused by the mutual friction. This possibility is particularly interesting since the instability 
may be intimately linked to the formation of vortex loops and superfluid turbulence (Andersson, Sidery & Comer 2007). In 
fact, the present analysis provides an important complement to our previous discussion of the turbulence problem. 

It is obviously necessary to ask whether this instability is likely to operate in neutron stars. For this to be the case, one 
would require the critical wavelength to be much smaller that (say) the size of the star. Otherwise, the plane-wave analysis 
does not apply. Prom our previous discussion (Andersson, Sidery & Comer 2007) we know that 

where ft ~ 2 x 10^'' cm'^/s and log(6/ao) ~ 20. From this we see that the critical wavelength for which the instability first 
appears is 



^^"250(^^^J cm-, (131) 
corresponding to a wavelength 

If we compare this to the typical intervortex spacing 



-1/2 



6«3.4xlO-^(^^^^J cm, (133) 

wo conclude that one may well expect modes with a wide range of wavelengths to be unstable in a typical neutron star. This 
is an interesting possibility, and it would be exciting to consider various scenarios where the instability may operate. 



7 TURBULENT MUTUAL FRICTION 

The presence of a dynamical instability in the vortex array will lead to oscillations in the vortices, triggering reconnections 
and the formation of vortex loops with a range of different sizes (Andersson, Sidery & Comer 2007). This behaviour is very 
similar to the standard cascade seen in normal fluid turbulence. If a turbulent tangle is present, then our analysis is no longer 
valid. After all, the form we arc using for the mutual friction force is based on the assumption that the vortex array is (locally) 
straight. One of the outstanding issues in superfluid Helium research concerns the nature of the force in the turbulent case. 
While some sort of consensus has been reached in the case of isotropic turbulence, problems with both relative flow and 
rotation are still far from understood. Yet this is the problem that we need to solve in order to model neutron stars. Our 
system is rotating, and if it becomes turbulent then any tangle that develops should be polarised. 

In absence of a clear strategy for developing a detailed model for the mutual friction force in the case of polarised 
turbulence, we have previously proposed a phenomenological prescription (Andersson, Sidery & Comer 2007). It is interesting 



Waves and instabilities in dissipative rotating superfluid neutron stars 15 



to apply this decription to the plane- wave problem we are currently investigating, since this may lead to a better understanding 
of the effect that turbulence may have on the vortex instability. 

We take as our starting point the mutual friction force posited in Eq. (78) of Andcrsson, Sidory & Comer (2007). In 
essence, this means that we add a term accounting for the presence of a turbulent tangle to (19) which now represents to 
polarisation of the vortex array. Thus we have 

= ... + ^nBwr , (134) 
where the . . . represents the straight vortex term from before. We have defined the total vortex length per unit volume. 



oi V aiPi W j.i/2 



4 V /3 



Lr , (135) 



where W = |WnpKi| = Vn, Lr = Uy and the constant parameters are such that (the arguments for this are given by Andersson, 
Sidery & Comer (2007)) 

with X of order unity. If we work out the perturbation of this new contribution to the overall force in the particular case when 
the wave vector is aligned with the background rotation, then wo find that SLr^ = 0. Moreover, when the imposed relative 
flow is also aligned with the rotation, we have dW = 0. This makes it very easy to account for this new force contribution. 
Under the conditions assumed in the previous section (protons clamped etcetera) , we simply get 

Sfr' = -..-^^Bv: . (137) 

When this term is added to the right-hand side of (122), and the dispersion relation is worked out as before, we find the wave 
solutions 

w = ±(2n + z/fe^) + iB{2a + vkl T V^kz) + ^^kB . (138) 

This suggests that turbulence always damps the inertial waves, as one might have expected. For the parameters given in (136) 
we see that 

This shows that the new damping term in (138) is very small in the neutron star case, when S <C 1. If we nevertheless include 
this contribution, and work out the critical velocities we find two roots. Assuming that kz is suitably large we retain (125) 
as the velocity at which the instability sots in. In addition, we find a second critical flow, beyond which the system is stable. 
The corresponding critical velocity is approximately given by 

K « ^Kfe . (140) 

For typical parameters, this velocity would be vastly greater than the critical velocity at which the instability sets in. In fact, 
it may well be the case that one can not reach such large relative flows in a realistic neutron star. Nevertheless, the result is 
conceptually interesting. One should also keep in mind that B is of order unity in superfluid Helium (Donnelly 1991), so this 
upper cut-off for the vortex instability may not be out of reach in that context. 



8 BRIEF SUMMARY 

We have analysed the wave propagation in a rotating superfluid neutron star cores, taking into account the standard mutual 

friction force. Our plane- wave analysis has added to previous discussions of this problem in a number of important ways. First 
of all, for models where the two background fluids co-rotate, we have clarified the role of chemical coupling and entrainment 
on both sound and inertial waves. Secondly, we have considered the mutual friction damping, demonstrating the well-known 
fact that sound waves propagating along a vortex array are undamped. Wo have also shown that the same is not true for 
inertial waves, which are damped by the mutual friction regardless of the propagation direction. We have accounted for the 
relatively small contribution of the vortex tension, which arises due to local vortex curvature. Focussing on purely transverse 
inertial waves, we derived the correction that the tension induces in the wave frequency. 

The most exciting result of our investigation concerns the presence of a dynamical instability associated with the inertial 
waves. The instability requires a linear relative flow in the background. We analysed the particular case when this flow is 
aligned with the rotation axis. This led to a demonstration that the mutual friction coupling induces an instability once the 
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relative velocity has reached a critical level. This instability is well-known from the analogous problem for superfluid Helium, 
and hence our result should not come as a great surprise. Nevertheless, the possibility that this instability may operate in 
neutron stars has only recently boon appreciated (Poralta ct al 2005, 2006; Andcrsson, Sidery & Comer 2007). We have argued 
that the instability belongs to the general class of two-stream instabilities. This interpretation is (we believe) new, and adds 
insight also into the Helium problem. 

If this instability operates in a neutron star, it is likely to lead to the formation of a vortex tangle and a state of superfluid 
turbulence. The impact of this on, for example, glitch recovery is not yet understood. Nevertheless, it is clear that much of 
our current "understanding" (which tends to be based on the assumption of a locally straight vortex array) may have to be 
revised. In view of this, the results we have presented here arc exciting. Having said that, it is clear that there are a number of 
difficult issues that need to be addressed if we really want to understand this problem. Our analysis was based on a number of 
simplifying assumptions, in particular we assumed that the proton fluid was clamped. It would be relevant to try to consider 
the general problem. One would certainly want to account for the entrainment, which will alter the critical velocity for the 
onset of the instability etcetera. It would also be relevant to try to quantify the damping (and possibly stabilising role) of 
shear viscosity, which should be important for short wavelength oscillations. We also need to consider various astrophysical 
scenarios for which the instability may be relevant. If it is the case that the key features required are a straight vortex array 
and some imposed relative flow, then the instability could be relevant in a number of situations. The most obvious possibilities 
would be i) neutron star free precession where the neutrons and protons essentially rotate with respect to diflferent axes (in 
the simplest model), ii) neutron star spin-down which (in a non-magnetic star) is faciliatcd via a viscous Ekman layer at the 
base of the crust inducing a global flow in the charged component, and iii) global mode oscillations, where the length scale of 
the mode is vastly larger than the typical lengthscale of the instability. These are all exciting problems, well worthy of further 
consideration. 
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APPENDIX A: THE VORTEX TENSION 

In this Appendix we provide the argument that leads to the form for the neutron vortex tension used in the main body of 
the paper. The calculation is based on the intuitive reasoning of, in particular. Hall (1958). It is important in the sense that 
it demonstrates how the entrainment parameters enter in the vortex tension. 

The starting point is the conservation of vorticity. Defining the macroscopic vorticity as 

i i 1 i 7 fc T-7 n / A 1 \ 

UJ —UyK — e VjPk , (Al) 

where the neutron momentum is Pn = Tin [v^ + £n(Vp — v^)] , we have 

/ UvKidV + / UvKiV^^dSj = . (A2) 
Jv Js 



D_ 
Dt 



IV Js 

Here it is assumed that the vortices move collectively with velocity v}^ . Use the divergence theorem to see that we must have 
dtoJi + VjicJivi) =0 . (A3) 
Now note that 

V.o;^ = , (A4) 
uj'Vjvi = . (A5) 

The first statement is trivial given the definition of the vorticity. The second should be true provided that there is no motion 
along the vortices themselves. This way the above conservation law can be recast as 

dtoji + (eijfee'°'"'w(V^) = . (A6) 
This leads to 

ezjfcV^ jStPn - e'°''"i^femnoV>n j = , (A7) 
which then requires that 
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ftp" - e^noV"K = V'=* , (A8) 

with "if some scalar potential. 

Let us now, for simplicity, assume that the only force that acts on the vortex is the Magnus force. Then we must have 

^yf+y^'i-- . (A9) 

The first term represents the smooth irrotational flow past the vortex, due to for instance the presence of all other vortices. 

The second term represents the self-induced flow that arises when the vortex is curved [see Appendix of Andersson, Sidery & 
Comer (2007)]. As wo have shown elsewhere, this term can be written 

ind.n f * 7 ^ Z t-t ^ k ~ ^ j ^ Zt-t ^ k / \ t r\\ 

V, =- Vtijkli-' K\lK =VtijkK'K\lK . (AlO) 

In order to use this vortex velocity in the equation of motion (A8), we note that 

klm I n . ind,n\ klm n , — lr-7 ^k /A11\ 

e \Vl +Vi ' \ Km = £ VlKm+ VlK . (All) 

Then we need 

klm ( n I ind,n\ nr-rfc I nr-rl k , — Ir-j ^k /A 1 o^ 

mnUye \Vi +Vi j Km = ViV Pn-ViV Pn + mnUyVK VlK . (A12) 

Use this in (A8) to get 

dtpt + vfVp^ - vfV''pl = V*'* + mnn.i^KVi/t'' . (A13) 
Finally use the definition of the momentum to get 

{dt + viVj)ivi + e„W*n) + EnW^nV^u" = V'x + UvVK^Vjfi' , (A14) 

where we recall that the velocity difference is Wpn = — This is the equation of motion for the superfluid neutrons, with 

the contribution from the vortex tension accounted for. The scalar potential x can easily be interpreted as the sum of the 
chemical and gravitational potentials to arrive at the standard form for this term. This way we arrive at Eq. (95) in the main 
text. 



APPENDIX B: THE VORTEX INSTABILITY IN A MORE GENERAL CONTEXT 

In this Appendix we provide a slightly more general derivation of the vortex two-stream instability that was dicussed in 
Section 6. While we still assume that the proton fluid is clamped, and neglect entrainment, we initially relax the assumption 
that the wave vector is aligned with the rotation axis. We also do not assume that the waves axe purely transverse. The results 
obtained in Section 6 follow in the appropriate limits, and the more complicated calculation that we outline here shows how 
the instability threshold can be derived under less constrained conditions. 

In the general case, the plane wave equation for the neutron fluid can be written, cf. (19), 

iCbVi + ifci/in + 2eijk^^v^ = (5/f , (Bl) 

where ui = u) + V^kj and we have used 5p,n = /Xne*^"*"'"'°^^^\ The force 5/° is made up of three contributions. The first is the 
mutual friction for a straight vortex array, and it leads to 

= Beijke'^'^ [ifliV^ie'^'kpV^ - Q.' e^^M^ k" vl) + ei^.k^ vlV^^, + 2^^^,^] ■ (B2) 

Next we have the contribution from the self-induced flow, which accounts for the vortex curvature. As long as the background 
flow is uniform, we can perturb (120) to get the contribution 

dfr"^ = vB{Cl^h)tijk(l' [e*"'"A;„i'" - n'=ep<,.n''fe<'t;^] . (B3) 
Finally, we have the vortex tension which is given by (102), i.e. 

= -u{h'k,) [eumk'v^: - hM^,(l'-k''vl)] . (B4) 
Putting all the pieces together and rearranging, the perturbed momentum equation can be written 
|w - iB[2n + i^{kjCl'f]jvl + [/in + iiyB{Cl'vf){Cl'^ki)]^ h - 2mjk^'vl 
= - BiV^Q")Clj - i:^{k„n")]e''''kkV^ - 2iBfl{Cl' v'})j Cli 

-B{ejimn'k'v:^)Vr + mVj'Cl^) + ii/(A;,n^)leii^fe'C ■ (B5) 
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Here the chemical potential perturbation is (in the clamped case we are also assuming that the proton density variation 
vanishes) given by 

P. = ^Pn . (B6) 

Since the continuity equation gives 

iuipn + ipn{kjvi) = , (B7) 
we get, using the standard definition of the sound speed from (35), 

Mn = -^(M^). (B8) 
The trick now is to form different scalar equations from (B5). By taking the scalar product with A;' we get 



2\ 



(kiVn) + ik,B(2Q. + vk'^)(Ct>v'^) 



UJ - iB{2n + vkl) -^{k^ - kl) 

U3 



+ [2in + ^^.(K'j^") + ivkl - B(fcnKr)] W = Bk,{t'''W]'kkvY) . (B9) 

In writing down this expression we have decomposed the wave vector into a piece along the rotation axis and a piece orthogonal 
to it, i.e., we are using 

k' = ki+ki , where k^Qj = . (BIO) 

We have also defined the scalar quantity W — eijk^^k-' v^. 

If we take the scalar product of (B5) with Cl wc find another scalar relation; 

i:j{0.'vi) - ^{kjvi) + B{Vj'Cl')W = B{e''''Vj'kkvf) . (Bll) 

UJ 

Prom the combination (B9)— fc^ x (Bll) we then get 

(kjvi) + [iBfe^(2n + i^k^) - k^uj] {Clivl) 

-\-[i{2Q. + vkl) + B{k,Vi)]W = Q . (B12) 
By taking the cross product between fc* and (Bo) wo get another useful relation. After some work it can be written 

[ik.{2a + vk") + Bk''{V^Ci')\v^ - [B{V^Cl') + ivk,]{hvi)ki 
= 2i{kjvi)Q.i - [w - iB{2Q. + iykl)]eimnk"'v'^ 
+ [[BVJ' - B{V^Q")nj - wk^Qjy''^kkvf - 2iBn{n'vj)^ eimpk'^CV 

-BWeipXV^ ■ (B13) 
Taking the scalar product of (B13) with f2' we arrive at 

[ikz{2n + vk^) + Bk'^{v^QP)]{Clivi) - [i{2Q. + vkl) + Bkz{Q.jVi)]{kivl) 

+[uj - iB{2n + vkl) + Bn]W = . (B14) 

Here we have defined yet another scalar Q, = eij^^vyv^ . 

We now have three equations, (B9), (B12) and (B14), for four unknown scalar quantities. To solve the general problem 
we need another relation. Although this relation can be obtained in a few steps by following the above strategy, we choose 
not to write it down here. Instead, we note that if we were to align with then^ 

= , (B15) 
eijkViiev'l = V^W . (B16) 

Then (Bll) simplifies to 

Cb{Cl,vi:)-^{k,vi)=Q . (B17) 

UJ 



^ We did not assume alignment from the beginning since we wanted to outline how the general problem would be solved. 
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We can use this relation to get an expression for kjvi. Using the result in (B12) and (B14) we only have two equations to 
solve. They are 

„2 



- iB(20 + lyki) 



+ kz[iB{2Q. + vk^ 



+[i(2n + vkl) + Bk^Vr,]W = 



and 



kzC^ 



[i(2n + vkl) + Bfe^Vn] + ikz{2Q. + vk^) + Z3fcV„ \ (fiiV, 



-]-[oj-iB{2Q. + vkl)\W = Q 
From this we easily obtain the dispersion relation 



(B18) 



(B19) 



- iB{2V. + vkl 



i{k' 



kl) 



+ kz\iB{2Q. + vk^) - w] - iB{2Q. + vkl)] 



-[{20. + vkl) - iBk.V^] - fez(2r2 + i/r) +iZ3fcV„ \ [{20. ^ vkl) - iBk.Vu] = 



(B20) 

It is, of course, not easy to write down the general solutions to this dispersion relation. But we can learn a lot from it if 
we make some further simplifications. To discuss these examples we first note that kz = k cos ^, where is the angle between 
the wave vector and the rotation axis. It is then straightforward to verify that we retain the solution from section 5 in the case 
when 9 = B = Q. We obviously also get the neutron sound waves. If we focus our attention on the possible vortex instability, 
then it would be natural to first relax the assumption that vanishes. Doing this, but still leaving B = Q (and in addition 
assuming slow rotation and weak tension) we find the leading order wave solutions; 



-2 

UJ = 



±kcn [l + {O sin^ e/k'^cl){20 + vk"^ cos^ 6)] , sound waves , 
±cos6l(2n + vk'^)^^^{20 + vk"^ cos^ , inertial waves . 



(B21) 



If we now linearise the dispersion relation in B, and assume that the modes take the form uj = uiq + BSQ, then we find that 
the mutual friction induced frequency correction follows from 



5oj = i 



Oo — kVn cos 



n sin^ 9 



sound , 



(B22) 



0(1 +cos2 e)+vk'^ cos'^e 



inertial 



Recalling that the waves are unstable if the imaginary part is negative, we see that (assuming that k cos 9 >0) the solutions 
for which wo < are always stable. In contrast, the uto > solutions become unstable at the critical velocity 



Vc 



Wo 



(B23) 



kcos9 

As one might have guessed, the onset of the instability depends on the projection of the wave vector along the relative fiow. 
For the sound waves we thus find that the critical flow is 

sound ^ ^n 



cos^ 



(B24) 



Since the superfluidity is likely broken before the wave propagation reaches the speed of sound, this indicates that these modes 
are always stable in a real system. Again the inertial waves are different. We flnd that 



-{20 + vk 
k 



2x1/2 



{20 + vk^ cos^ 6*) 



1/2 



(B25) 



According to this criterion, the instability actually sets in at a lower relative velocity when the wave vector is not aligned 
with the rotation axis. Of course, in reality one may expect the tension term to be small compared to the rotation term. Then 
the difference between the above result and the aligned case discussed in Section 6 may only be signiflcant at extremely short 
wavelengths. 
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